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Abstract. We present a comprehensive investigation of 
the radio polarization properties within the theory of 
natural electromagnetic wave modes in pulsar magneto- 
spheres. Taking into account the curvature of the field 
lines, aberration effects and magnetic sweep-back we use 
the relativistic dielectric tensor in the low-density cold 
plasma approximation and derive the following polariza- 
tion characteristics, which are in full agreement with the 
observational findings. Specifically, we demonstrate that 
1. The degree of linear polarization decreases with in- 
creasing frequency. 2. The degree of circular polarization 
increases with increasing frequency. 3. At high frequencies 
(> a few GHz) the degree of linear polarization is corre- 
lated to the spin down luminosity E. 4. At high frequen- 
cies long-period pulsars exhibit weaker linear polarization 
than their short-period counterparts. 5. The difference be- 
tween the refractive indices of the two natural wave modes 
decreases with increasing frequency which possibly results 
into a depolarization via superposition. 



Key words: Plasmas; Polarization; Radiative transfer; 
Waves; Pulsars: general 

1. Introduction 

The polarization properties of pulsar radio emission are 
highly diverse. Nearly all possible polarization states have 
been detected, from pulsar to pulsar, from frequency 
to frequency, from pulse to pulse and even within one 
pulse profile. Nevertheless, some general statements can 
be made from an observational point-of-view, which must 
be explained by any plausible theoretical model: 

— Pulsar radio emission is usually highly linearly pola- 
rized at low frequencies (in the following "low fre- 
quencies" defines frequencies up to about one GHz, 



"high frequencies" corresponds to everything well 
above one GHz). Towards higher frequencies the ra- 
diation tends to depolarize (e.g. Manchester 1971 



Xilouris et al. 1996) 



— The observed polarization position angle (hereafter 
PPA) is a very stable feature which usually does not 
depend strongly on frequency. Apart from the occur- 
rence of sudden ~ 90° jumps, the so-called orthogonal 
polarization modes (hereafter OPM) , the PPA usually 
follows a characteristic S-shaped curve. This curve can 
well be modelled by the rotating vector model (here- 
after RVM) as proposed by Radhakrishnan & Cooke 
(1969). It is therefore thought to reflect the geometry 
of the radiating field lines. 

— The sum of many (typically at least a few hundred) 
single pulses usually converges towards a characteristic 
average profile. In contrast, the individual pulses show 
a high variability. One well known single pulse fea- 
ture is the quasi-periodic microsecond- structure (e.g. 
Cordes &: Hankins 1977 , Lange et al. 1998 ) observed 
over a broad frequency interval with no obvious fre- 
quency dependence of the quasi-periodicity. It is usu- 
ally highly linearly and circularly polarized with a 
roughly constant PPA and often OPMs at the edges. 
This suggests that, during one micro-pulse, only one 
of the orthogonal polarization modes is observed. An- 
other feature are the subpulses, sometimes drifting 
in pulse phase from one pulse to the next one. The 
PPA usually shows a swing during one subpulse, with 
the swing being stable relative to the subpulse. So 
when the subpulse drifts, the swing drifts as well, 
causing a depolarization in the integrated profile (e.g. 
Manchester et al. 1975). 



Send offprint requests to: A. von Hoensbroech (avh@mpifr- 
bonn.mpg.de) 



In average profiles and in individual pulses OPMs can 
be observed (e.g. Stincbring et al. 1984, Gangadhara 
1997). The magnitude of these jumps is usually 90° but 
can also be less sometimes. Sudden orthogonal jumps 
in the PPA can be found in nearly all pulse phases 
and even at the same pulse phase in successive sin- 
gle pulses. The OPMs often coincide with changes in 
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the handedness of the circular polarization. Their ex- 
istence therefore suggests that the radiation preferen- 
tially takes two orthogonal elliptical states. 
Some pulsars show a high degree of circular polar- 
ization, preferentially near the centroid of the pro- 
file. Sometimes it also sometimes shows a sense rever- 



sal near the centroid ( Radhakrishnan fc Rankin 1990 ) . 
Recent observations show that a few pulsars have a 
strong increase in their degree of circular polariza- 
tion towards high frequencies together with a decrease 
in linear polarization, thus suggesting that a process 
might be active which transforms linear into circular 
polarization ( von Hoensbroech ct al. 1998 ). 
At low frequencies there seems to be no signifi- 
cant relation between the average degree of polariza- 
tion and any other pulsar parameter. At higher fre- 
quencies, however, correlations have been found be- 
tween the degree of polarization and the spin down 
luminosity E respectively the surface acceleration 
potential (Morris et al. 1981, Xilouris ct al. 1995| and 
von Hoensbroech et al. 1998). Pulsars with a high E 



(respectively a short period and a large period deriva- 
tive) show a higher average polarization at high fre- 
quencies than those with a lower E. This suggests 
that for these pulsars a possible depolarization pro- 
cess becomes important at higher frequencies than 
for the other pulsars. Closely connected to this is an 
anti-correlation between the degree of polarization at 
high frequencies and the period. Only rapidly rotating 
pulsars seem to have highly polarized radio emission 
( von Hoensbroech ct al. 1998 ). 

Several authors estimated the emission altitudes where 
the observed radiation is produced and studied if there 
is a frequency dependence of this altitude (Cordes 
1978, Rankin 1990, Thorsett 1991, Blaskiewicz et al. 
1991, Gil 1991, Phillips 1992, Gil & Kijak 1997, von 
Hoensbroech & Xilouris 1997, Kramer et al. 1997). 
Various independent methods were used and all au- 
thors came to the conclusion that the radio emission 
originates from a region close to the pulsar surface at 
a few percent of the light cylinder radius (at least for 
"normal" - non millisecond - pulsars). The frequency 
dependence of the emission height is thought to be 
small, if it exists at all. 



All these observational facts have to be explained 
within a general model for the polarization of pulsar radio 
emission. Basically, there are two aspects which must be 
considered. One is the polarization characteristic of the 
emission process itself, and the other one is the influence 
of the magnetosphere on the propagating radiation. In this 
paper we concentrate on the possible role of propagation 
for the polarization because this can be considered inde- 
pendently of the possible radiation mechanisms. 

The magnetosphere of a pulsar is filled with a plasma 
which is streaming relativistically outwards along the open 



field lines. A number of articles has already been published 
on propagational effects in this plasma. 

Cocke & Pacholczyk (1976) considered the effect of 
Faraday pulsation for quasi-transverse propagation. Fara- 
day pulsation is the general case of Faraday rotation. It 
occurs when a polarized beam decays into two elliptical 
orthogonal modes which propagate with different phase 
velocities. The PPA of the resulting beam gets rotated 
and the polarization changes between linear and circular. 

In a series of papers Melrose & Stoneham (1977), Mel- 
rose (1979) and Allen & Melrose (1982) derived an ap- 
proximation for the dispersion tensor and calculated the 
properties of the natural wave modes in the plasma (see 
also Lyutikov 1998). They assumed that the polarization 
of the propagating electromagnetic wave follows the shape 
of the natural wave modes up to a certain radius of lim- 
iting polarization (hereafter i?Lp)- The radiation then es- 
capes with the polarization properties of the natural wave 
modes at -Rlp- Barnard (1986) places -Rlp at the place 
of cyclotron resonance. For normal pulsars this resonance 
occurs at a few 10% of the light cylinder radius, but for 
short period pulsars (P < 0.06s) this resonance is outside 

Cheng & Ruderman (1979) considered adiabatic walking 
to be responsible for the variable polarization in sub- 
pulses. The polarization properties were derived for dif- 
ferent emission mechanisms and plasma compositions. 

Harding & Tademaru (1979, 1980, 1981) presented numer- 
ical calculations on the propagation of linearly polarized 
pulses through a shearing plasma. They show that this 
can account for micro-structure, circular polarization and 
rotation of the PPA. 

Barnard & Arons (1986) and Arons & Barnard (1986) 
have calculated the dispersion relations for the X- and the 
O-mode in an ultra-relativistic, one-dimensional plasma 
for distribution functions in e~ and e + (the O-mode has 
its electric field vector in the plane of curvature of the 
magnetic field, the X-mode perpendicular to it). The X- 
mode is purely linearly polarized and propagates easily 
through the magnetosphere whereas the O-mode can have 
some circular contributions and follows the bending of 
the field lines. This causes a separation of the OPMs 



flMcKinnon 1997p . 

Calculations concerning the propagation characteristics of 
the various natural modes have also been presented by Be- 
skin et al. (1993). Using the limit of infinite magnetic field 
strength they conclude that only the electromagnetic X- 
mode can propagate freely through the magnetosphere. 
The angular dependence of the properties of the modes 
and the influence of the cyclotron resonance are derived. 
This resonance becomes only significant, if the secondary 
particle production rate is high (> 10 4 per primary par- 
ticle) thus leading to a much higher plasma density. The 
different mode properties also account for the circular po- 
larization often observed in core components. 



von Hoensbroech et al.: Natural Polarization Modes in Radio Pulsar Magnetospheres 



3 



In this paper we extended the ideas proposed by Mel- 
rose & Stoneham(1977), Melrose (1979) and Allen & Mel- 
rose (1982). We derive the properties of the natural wave 
modes throughout the magnetosphere with special respect 
to the angle between the propagating wave and the mag- 
netic field. We make predictions for the qualitative de- 
pendence of the polarimetric properties with frequency 
and various pulsar parameters. These predictions are then 
compared with the observations. 

2. Method 

2.1. Global picture 

For our calculations we use the following global picture 
for the pulsar radio emission (e.g. Manchester & Taylor 
1977): The neutron star (hereafter NS) has a strong dipo- 
lar magnetic field which produces a corotating magneto- 
sphere filled with charged particles which are drawn from 
the NS surface. It is convenient to define the light cylinder 
radius i?LC of the magnetosphere at which the corotation 
velocity is equal to the speed of light 



-Rlc 



cP 
2^ ' 



(1) 



where c is the speed of light and P is the rotation period 
of the NS. Those field lines which do not close within i?Lc 
define the polar cap. The angular size of the polar cap is 
therefore given by 



' Rss 



(2) 



with i?NS being the neutron star radius (i?NS = 10 km, lof 
stands for last open field line). The particles get acceler- 
ated above the polar cap and flow relativistically outwards 
along the open field lines. The motion is one dimensional 
(all particles are in their lowest Landau-levels). Only one 
type of particles is considered (here electrons). Some pair 
producing processes might be active, thus adding another 
type of particles but we assume that there will always be 
one type dominating. The problems of strong pair produc- 
tion follows from considerations by Kunzl et al. which will 
be presented in a forthcoming paper. 

According to the model of Goldreich & Julian (1969) 
we use the following density: 



47re B Q 
e ' P 



(3) 



with e being the elementary charge and eo — 8.85 • 
10~ 12 AsV _1 m the electric field constant. The surface 
magnetic field is defined as 



B n = 



3c 3 hqI 



327T 3 i?6 s sin 2 a 



■ PP . 



(4) 



Here I = 10 38 kg • m 2 is the moment of inertia of the NS 
and a is the angle between the axis of rotation and the 
magnetic field axis. P is the temporal derivative of the 
period, (j,q = 4tt • 10 _7 VsA _1 m~ 1 the magnetic field con- 
stant. As the field is dipolar, the magnetic field strength 
and the local plasma density decrease with the cube of the 
distance from the NS. 

At a few percent of i?LC the plasma produces radio 
emission. Since we want to make our calculations inde- 
pendent of the actual emission mechanism, we only as- 
sume that the particles initially have a certain Lorentz 
factor 7bcam, just before they radiate. The angle of emis- 
sion is approximately 1 / 7beam with respect to the field line 
along which they stream. After having radiated the parti- 
cles flow further outwards and form the background plasma 
(bg) through which the radio waves have to propagate. 
Naturally the background plasma has a lower Lorentz fac- 
tor than the radiating beam, i.e. 7b g <C Tbeam- Since the 
radio emission mechanism must be coherent and therefore 
very efficient, we assume that the energy spread of the 
background plasma is negligible above the emission re- 
gion, thus the plasma is cold. In this case the cold plasma 
approximation applies and the energy distribution func- 
tion is 



/(l) - <K7-7bg)- 



(5) 



Further, we assume that we observe only one natural 
wave mode at a time. The polarization properties of the 
wave correspond to the local form of the natural wave 
modes up to -Rlp- Above this radius, the interaction be- 
tween the propagating wave and the plasma becomes so 
weak that the escaping wave preserves the state of polar- 
ization of that radius to infinity (see e.g. Melrose 1979, 
Barnard 1986 and Beskin et al. 1993). The maximum 
value for i?LP is i?LC but for most pulsars i?LP is probably 
reached at a closer distance. The minimum value for i?Lp 
is the emission height itself. 

The assumption that only one polarization mode is 
observed at a time might not always be true, but the ex- 
istence of highly polarized OPMs indicates that this is at 
least sometimes true. If both modes are observed simulta- 
neously the radiation gets depolarized. The results of the 
following discussion can therefore be regarded as upper 
limits for the observed degrees of polarization. 

2.2. Geometrical Aspects 

It is very important to analyze carefully the angle of the 
propagating wave relative to the magnetic field. This can 
be understood by an inspection of the extreme cases for 
parallel and transverse propagation. For an electromag- 
netic wave propagating parallel to a magnetic field, the 
natural wave modes are circular (R- and L-mode) , whereas 
for transverse propagation they are linear (X- and O- 
mode) (see e.g. Krall & Trivelpiece 1986 or any standard 
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plasma physics book). For oblique propagation intermedi- 
ate states appear (see Fig. ||). It is our aim to calculate 
these states as a function of the various parameters in- 
volved. 

To derive the angle between the propagating wave and 
the magnetic field lines the following contributions have to 
be considered: The angle which originates from the dipolar 
structure of the field, the angle associated with the aber- 
ration and the degree of magnetic sweep back. Aberration 
is the apparent forward bending of the field lines through 
the rotation of the pulsar ( Phillips 1992) ) . The magnetic 
sweep back is caused by the magnetic dipole radiation. 
The magnetic field lines are bent in the direction opposite 
to the pulsar rotation by the electromagnetic torque of the 
neutron star Qhitov 1983| ) . 

We do not consider the relativistic beaming angle 
(^beaming = l/7bcam) outside the emitting region itself be- 
cause this angle is symmetrical with respect to the direc- 
tion of emission and will be averaged out statistically. 

2.2.1. Dipolar contribution 



field line 




x(x, e n ) 



@(x,e ) 



neutron star 



Fig. 1. Geometry to derive the tangential inclination x{ x ^o) 
of any dipolar field line with the footpoint-angle Oq at the radial 
distance x. 



In a polar coordinate system, the quotient sin 2 9/r is con- 
stant for a dipolar field line. Having the NS in the centre, 
and defining each field line by its footpoint (point of inter- 
section) at the NS-surface (i?NSi$o)i a parametric equation 
describing the field line is given by: 



= constant 



sin 2 9 
Rns 



(6) 



Some straightforward calculations lead to an expression 
for the tangential angle x( x , 0o) of the 6*o-field line: 



x{x, 9q) = 8(x, #0) + arctan - tan 9(x, #0) 



1 



Here x — 7'/i?NS and Q(x,9q) = arcsin( v / x • sin#o) is the 
radial angle of the ^o-field line at a distance x (see also 
Lync & Graham-Smith 1990). The geometry is displayed 
in Fig. [l]. This can be used to calculate the polar angle 9s 
for a given distance x at which a ray propagates, that was 
emitted tangentially from the 9o-&eld line at the emission 
height x cm : 



9s(x,6o) = x(x om , #0) - arcsin 



X{ x em, 9q) + 6(x em , #o)) 



■ sin (jr- 



(8) 



Using this point (x,9s), the footpoint-angle 6*so of the cor- 
responding field line can be calculated using Eq. (^|). 



tangentially emitted beam 



Am 




emitting 
field line 



-x(x,e so (x,e )) 

«< v em' y 



neutron star 

Fig. 2. Geometry to derive the intersecting angle ipi(x, #0) be- 
tween a beam which was emitted tangentially from the 9o-&e\d 
line at the emission height x cnl and a field line in a radial dis- 
tance x. 



The intersecting angle between a ray emitted tangen- 
tially from the 9q field line in the emission height x cnl with 
the field line at the distance x can then be found as the 
difference of the tangential angles 

6> ) = x(x, 6so{x, Oo)) ~ x(^om, 0o ) , (9) 
as is shown in Fig. |^. 

2.2.2. Other contributions 

The rotation of the magnetosphere causes additional ef- 
fects which have to be considered for the intersecting an- 
gle. The radiation is bent forward into the direction of 
rotation by aberration. Following Phillips (1992), the mag- 
nitude of this effect is 



(7) a b(^) = arctan ( 2-7T 



xi?NS sin a 
~~cP 



(10) 
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At larger distances from the NS, corotation of the 
plasma cannot be maintained, and magnetic sweep-back 
becomes important. The calculation of this effect is rather 
complex, but an approximation was given by Shitov 
(1983): 



9 msh (x) ~ 1.2 



/ 2ttx -Bns \ ' 



sin 2 a 



(11) 



The direction of this effect is backwards relative to the 
rotation. To calculate the net-effect of these two contri- 
butions, one has to take their difference at the distance x 
and subtract their value at the emission height x cm . It is 
therefore given by 



4>2{x) = 6 ah (x) - 6> ab (Xom) 



+ msh (x cm ) . (12) 



Angular contributions by effects of general relativity are 
negligible at distances larger than a few km from the NS 
surface (e.g. Gonthier & Harding 1994). Since the contri- 
bution will always be much smaller than the accuracy of 
our calculations, which is restricted to the beaming angle 

Of l/7bcam- 

2.2.3. Resultant intersection angle 



into account. This can be defined as the angle between the 
plane of curvature of the emitting field line and the field 
line in the middle of the profile cV The tube of emitting 
field lines then rotates into the direction of the field line 
with 0o = 90°. 

The resultant angle of a beam, radiated from the 9q- 
field line into the direction cf>o at the emission height x em , 
which intersects with a field line in a distance x from the 
NS can then be given as: 



V'sOMo) = [(V>i(x,6> )cos</>o) 2 + 

+ (Mx)- ^Mo)sin<ft)) 2 ]' 



(13) 



This equation is only valid if one regards the inter- 
secting angle at a sufficient distance from the emitting 
altitude, where tps reaches a finite value. At the emission 
region itself, ips would be zero. As the radiation is beamed 
into a forward-cone of l/7bcam, and the intersected field 
lines are all nearly parallel in the emission region, we use 
this beaming-angle ip s = 1/Tbeam for x ~ x cm . 

In Fig. |i[ the intersection angle between a tangentially 
emitted beam and a field line at a given distance is dis- 
played from the emission height to i?LC for a canonical 
pulsar. 



Emission cone, seen from the top\ ,■ r . 

r \ line ot sight 




Contribution through Contribution through/ 
dipolar structure aberration and magnetic sweepback 



2.2.4. The canonical pulsar 

As we require a number of intrinsic pulsar parameters 
for our calculations, we use the mean values from the 
known pulsar sample. These parameters are P — 0.6 s 
and P = 10~ 146 . For the magnetic inclination we as- 
sume a = 40°. The regarded field line is at cf>o — 40° and 
#o = 0iof/l-5 (8\ f is the footpoint angle of the last open 
field line, see Eq. ||). As a typical emission height we choose 
2% of i?LCj which corresponds to slightly more than 50 
pulsar radii for the canonical pulsar. For the Lorentz factor 
of the background plasma we use be 7b g = 1.7, following 
Kunzl (1997). In those cases, where we derive properties 
at a given distance in the magnetosphere, we use a point 
at 20% of the light cylinder radius for i?LP (see also Sect. 
3.1). As a typical value for Tbeam we assume a value of 50 
(e.g. Rathnasree 1996| , Kunzl 1997). As stated earlier, the 
value of 7beam is only important very close to the emission 
region itself. 



Fig. 3. Geometry to derive the resultant intersection angle 
under consideration of the dipolar contribution and those of 
aberration and magnetic sweep-back, tpi points towards the 
direction of the emitting field line, # a b into the direction of 
rotation. 9 ms b points backwards relative to the direction of ro- 
tation. 

For the resultant intersection angle (see Fig. |||), the cur- 
vature direction of the emitting field line has to be taken 



2.3. Approximation for the Dielectric Tensor 

In order to derive the properties of the natural wave modes 
it is necessary to get an approximate idea of the dielectric 
characteristics of the plasma. Our approach follows the 
work of Melrose (1979), which is based on two assump- 
tions: 

— The plasma is in the low density limit, which is defined 
by the condition that the Alfven velocity is much larger 



(i 
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0.0 0.2 0.4 0.6 0.8 1.0 

x: Distance from NS (1 = light cylinder radius) 

Fig. 4. Resultant intersection angle for a canonical pulsar (see 
text for assumed parameters) . The angle is displayed from the 
emission height to the light cylinder radius. 

than the velocity of light (v\ 3> c). Combining the 
Alfven velocity 
B 2 



with the relativistic mass density p 
to the condition 
B 2 

-5 > 1 . 



7b g r 



(14) 
n leads 

(15) 



&■ ■ n Q ~fm c n 

This condition is well fulfilled throughout the pulsar 
magnetosphere. For our calculations this implies that 
the dielectric tensor is close to unity. 
— The background-plasma is assumed to be cold. This 
follows from the consideration that the coherent ra- 
diation mechanism can only be effective for particle 
energies above a certain threshold (see Eq. ||). The fi- 
nal particle energies will therefore peak at this critical 
value. 

Thus the long wavelength approximation applies, which 
can be used when the wavelength of the propagat- 
ing wave is much larger than the relevant lengths in 
the plasma. This approximation should be fulfilled 
throughout in the pulsar magnetosphere above the 
emitting region. 

Using the low density limit the dielectric tensor may 
be written as 



Ae, : 



(16) 



with Aei, <C 1. The dispersion equation is then found as 



An 4 - Bn 2 + C = 



(17) 



with A, B and C being functions of ey. In this approxi- 
mation Aeij does n ot depend on n. Hence two solutions 
can be given for n 2 ( Melrose fc Stoncham 1977 ): 



i 2 = l±Anl 



An 2 = 



+ b±^(a-b) 2 + 



(18) 
(19) 



Using the cold plasma approximation for electrons the di- 
electric tensor leads to 



^pe sin2 1 



n 2 



b = - 



w 2 (l -/3 cos 6>) 2 

io 2 c (l-/3cosfj) 2 
w 2 (1-/3 cos (9) 2 -fl 2 
(3) 



w 2 (f2 e /w)(cos( 



lo 2 {\ — (3 cost 



n 2 ' 



(20) 
(21) 
(22) 



(3 = v/c. The electron plasma frequency 07 pc and the elec- 
tron gyro frequency fi c arc defined as: 



■ipe 



a. 




7m 



(23) 
(24) 



n is the plasma density (n(x) = uqj/x 3 , see Eq. ||), e 
the elementary charge, m e the electron mass and eo the 
electric field constant. B is the magnetic field strength 
(B{x) = B /x 3 , see Eq. |). 

With these parameters, the axial ratio of the polariza- 
tion ellipses of the two natural wave modes can be calcu- 
lated: 



T+ 



2<? 



b Ty /(a-b) 2 + 4g 2 



(25) 



T = and T = ±oo correspond to linear polarizations. 
For elliptical polarizations T becomes finite. The rotating 
sense of a positive T is right hand with respect to the 
direction of propagation and left hand for negative T. 

The degree of circular polarization can then be calcu- 
lated through 



IIv± := sin(2 arctanT-t) 



(26) 



As a natural wave mode is fully polarized per definition, 
the linear polarization can be calculated as the quadratic 
complement: 



n 



L± 



1 



TT 2 
11 V± 



(27) 



3. Results 



The polarization of the natural wave modes depends 
mainly on the angle between the propagating electromag- 
netic wave and the magnetic field. If the wave propagates 
purely parallel to the field (k || B, k being the wave vector 
of the propagating beam), the modes obviously have to be 
circularly polarized since no linear direction perpendicular 
to k is preferred. When k is perpendicular to B, a direc- 
tion is preferred and the wave is fully linearly polarized. If 
the angle between k and B becomes larger than 90° then 
the handedness of the circular polarization gets reversed 
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Linear polarization 



Left hand 



Right hand 



Circular polarization 



45 90 135 180 

Angle between wave vector and magnetic field [deg] 

Fig. 5. Degree of polarization of one natural wave mode versus 
the intersecting angle between the propagating wave and the 
magnetic field. In this figure we show the principal polariza- 
tion behaviour of one natural mode of a wave propagating at 
an oblique angle to the magnetic field through a non relativis- 
ts (7b g = 1) cold electron plasma. The function was calculated 
for a frequency of 1 GHz and a magnetic field strength of 0.2 
T. For small angles between k and B the polarization is mainly 
circular, as no linear direction is preferred perpendicular to k. 
At an angle close to 90° the polarization is mainly linear, for 
angles larger than 90° it becomes increasingly circular again 
but with an reversed handedness (due to the opposite direc- 
tion of gyration). The transition between fully circular to fully 
linear depends heavily on the frequency and the magnetic field 
strength. For plasma Lorentz factors 7b g > 1, the angle gets 
Lorentz transformed - thus shifting the condition for transverse 
propagation to smaller angles. 



for each mode due to the opposite direction of the particle 
gyration. 

The transition (for an example see Fig. ^ between 
these two states depends strongly on a number of para- 
meters: 

— The magnetic field strength. For a strong magnetic 
field, the transition from circular to linear polarization 
occurs for very small angles between k and B. 

— The frequency. For lower frequencies, a higher degree 
of linear polarization is reached at smaller angles than 
for higher frequencies. 

— The plasma 7-factor. For a relativistic plasma, the in- 
tersecting angle between the propagating beam and 
the magnetic field is Lorentz-transformed, such that 
the angle increases. The transverse propagation (~ 90° 
in Fig. ||) is reached at a smaller angle. 



3.1. Qualitative polarization characteristics 

We apply these calculations to the propagation of an elec- 
tromagnetic wave through the pulsar magnetosphere using 
canonical pulsar parameters (Sect. 2.2.4 ) and taking into 
account all angular effects. 




0.2 0.4 0.6 O.f 

x: Distance from NS (1 = light cylinder) 



Fig. 6. Degree of linear polarization of one natural wave mode 
versus the distance from the NS (u p to t he light cylinder ra- 
dius) for the canonical pulsar (Sect. 2.2.4). 



Figure ^ illustrates the change of the natural wave 
mode polarization characteristics with distance from the 
pulsar for different frequencies. For clarity, we displayed 
only the linear polarization for one mode (the + mode, 
see Eq. 18). As the circular polarization is given by the 
quadratic complement (see Eq. ^7]) , a decrease of the linear 
polarization leads to an increase of the circular contribu- 
tion. The dominating parameter is the rapidly decreasing 
magnetic field. After the radiation, the kinetic energy of 
the plasma is relatively low. Following Kunzl (1997), we 
use a background Lorentz factor of 7b g = 1.7. This value 
would result from a free electron maser emission mecha- 
nism through comparing bunching- and loss times. Other 
radiation mechanisms lead to different values of 7b g but 
the qualitative manner does not change. The dependence 



on 7b g is discussed separately in Sect. 3.2.3 



The following qualitative behaviour can be observed. 
Close to the emitting region, the mode is fully linearly 
polarized for all frequencies due to the extremely high 
magnetic field in the emitting region (|jB| ~ 10 3 T). As the 
wave propagates outwards, the intersecting angle increases 
and the magnetic field decreases. Initially this just causes 
the linear polarization of the high frequency waves to de- 
crease, while the lower frequency waves maintain a high 
degree of linear polarization up to a large distance from 
the NS. Thus, if only one mode is observed at a time, for a 
given distance we expect a decrease of linear polarization 
with frequency and an increase of circular polarization. 

As stated earlier, it remains unclear up to which dis- 
tance from the NS the polarization of the propagating 
radiation is influenced by the changes of the natural wave 
modes (-Rlp)- This remains to be a problem to be solved 
by mode coupling theory but appears to be rather compli- 
cated as the plasma is most likely highly inhomogeneous 
in space and time. Therefore for the present discussion 
we choose an arbitrary distance for i?LP and derive the 
qualitative behaviour without loss of generality. 
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Fig. 7. Degree of polarization of one natural wave mode ver- 
sus frequency. The height above the NS was chosen to be 20% 
of the light cylinder without loss of generality. The qualitative 
behaviour is a decrease of linear and an increase of circu lar po - 
larization with frequency for the canonical pulsar (Sect. 2.2.4). 



Figure shows the frequency development of one nat- 
ural wave mode at an arbitrary distance of 20% of Rlc for 
a typical pulsar. The degree of linear polarization should 
decrease and the degree of circular should increase above 
about one GHz. 

3. 2. Dependence on various pulsar parameters 

We now want to investigate, how the change of polariza- 
tion characteristics with frequency depends on the varia- 
tion of other pulsar parameters. 

3.2.1. Spin-down luminosity 

As stated in Sect. |[ a correlation between the degree of 
polarization and the spin-down luminosity E has been ob- 
served at high frequencies. 



E = Att z 



P 



(28) 



I is the moment of inertia. In order to perform our cal- 
culations for different values of E, we keep all parameters 
fixed and change only the time derivative of the period 
P. As before, the period was set to 0.6 s and for P we 
used the values lO" 16 3 , lO" 14 6 and 10" 121 . This corre- 
sponds to E values of 9 • 10 23 W, 5 • 10 25 W and 1 • 10 28 
W respectively 

Figure || shows that for high E pulsars the degree of 
linear polarization decreases at a higher frequency than 
for low E pulsars. This is in agreement with the observa- 
tions, especially if one follows the argument of Sect. 2T 
and takes the degree of polarization as an upper limit. At 
low frequencies, all pulsars have a high degree of linear 
polarization as an upper limit. Therefore no correlation 
with E is expected. At higher frequencies - in contrary - 



Fig. 8. Degree of linear polarization of one natural wave mode 
versus frequency for different spin down luminosities E. We 
used the same canonical pulsar parameters as before except 
that we varied the values for P, keeping the period fixed at 
P = 0.6 s. The resulting values for E are written next to 
the curves. For pulsars with a higher E the decrease in linear 
polarization occurs at a higher frequency. 



Correlation coefficient = 0.1 



! 'I 



i > I 



H 1 I I 



26 

log (Edot [W]) 



Correlation coefficient = 0.8 



5 f' 1 I 
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Fig. 9. Measured degree of linear polarization versus the spin 
down luminosity E at 408 MHz (upper plot) from Gould & 
Lyne (1998) and 4.9 GHz (lower plot), taken with the 100 m 
Effclsberg radio telescope. The two plots consist of (nearly) the 
same set of pulsars. These observations correspond to the pre- 
dictions of Fig. [| At the higher frequency, pulsars with a low 
E have a low degree of linear polarization. At low frequencies 
no such correlation can be observed which again corresponds 
to Fig. ^. At low frequencies the degree of linear polarization 
is nearly independent of E. 
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the upper limit for low-E pulsars is low, whereas high- 
E pulsars can still exhibit highly linearly polarized radio 
emission. A correlation is therefore expected at those fre- 
quencies. 

3.2.2. Period 

The next parameter we vary is the period. In the list of 
observational constraints in Sect, [l] it was stated that only 
for short-period pulsars highly linearly polarized high fre- 
quency radio emission has been observed. We therefore 
keep all parameters fixed and vary only the period. The 
emission height was set to 2 % of the light cylinder radius 
(and therefore changes in absolute numbers as the light 
cylinder radius varies with the period). 

Figure [l^ shows that short-period pulsars can main- 
tain a higher degree of linear polarization up to a higher 
frequency than long-period pulsars. This is in excellent 
agreement with the observations. At low frequencies the 
linear polarization can be high for all periods, we there- 
fore would not expect an observable correlation. At high 
frequencies in contrast we therefore do expect an inverse 
correlation. 




Frequency [Hz] 



Fig. 10. Degree of linear polarization of one natural wave mode 
versus frequency for different pulsar periods. P was held fixed 
at 1CT 14 - 6 , the emission height at 2% and i? L p at 20 % of i? L c- 
Short-period pulsars can maintain a higher degree of polariza- 
tion up to a higher frequency than their long-period counter- 
parts. 



Correlation coefficient = -0.1 



3.2.3. Lorentz-factor of the background plasma 

The Lorentz-factor of the background plasma is important 
because the intersecting angle gets Lorentz transformed. 
If one increases the Lorentz factor at a given point, the 
Lorentz-transformed intersection angle first increases up 
to 90° in the plasma frame, resulting in fully linearly pola- 
rized modes of transverse propagation (see Fig. ||) . A fur- 
ther increase of 7b g leads to angles > 90°. The degree 
of linear polarization then decreases again and the circu- 
lar contribution changes its handedness. This behaviour 
can be seen in Fig. |l2|. One consequence is that if the 
Lorentz transformed intersecting angle between the prop- 
agating wave and the magnetic field is close to perpendicu- 
lar, small, local variations in the conditions can result into 
changes in the handedness of the circular polarization. As 
we have to assume that the plasma conditions are highly 
inhomogeneous in space and time, this effect could well ac- 
count for sudden changes in handedness observed in low 
degree circular polarizations. 

An inference which can be made is that the background 
Lorentz factor cannot be very large, e.g. larger than a few 
hundred. 



3.3. Refractive indices 

So far we only regarded the polarimetric properties of one 
natural wave mode. But, as mentioned in Sect. |], obser- 
vations suggest that both modes are present in the pulsar 
radio emission as the familiar OPMs. They can be ob- 
served both separated (leading to highly polarized parts 
in the pulse profile with a smooth variation of the PPA) 



■S 60 - I 



Period [s] 



Correlation coefficient = -0.6 



1.0 
Period [s] 



Fig. 11. Measured degree of linear polarization versus period 
at 408 MHz (upper plot) from Gould & Lyne (1998) and 4.9 
GHz (lower plot), taken with the 100 m Effelsberg radio tele- 
scope. The two plots consist of (nearly) the same set of pulsars. 
This observations correspond to the predictions of Fig. At 
the higher frequency, pulsars with a long period radiate with 
a low degree of linear polarization. At low frequencies no such 
correlation can be observed which again corresponds to Fig. [lo| 
At low frequencies the degree of linear polarization is nearly 
independent of the period. 
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Frequency [Hz] 

Fig. 12. Degree of linear polarization of one natural wave mode 
versus frequency for different background Lorentz-factors (the 
values are written next to the curves). As the intersecting an- 
gle between the propagating wave and the magnetic field is 
Lorentz transformed, the degree of linear polarization first in- 
creases with 7b g until the transformed angle reaches 90° and 
then decreases again for larger resulting angles. The difference 
is the handedness of the circular contribution for high and low 
7b g factors. 



O-mode 



10" i1 



Frequenz [Hz] 

Fig. 13. Deviations of the two squared indices of refraction 
of the natural wave modes from unity (n 2 — 1 + An±). Note 
that An 2 , is negative. An+ (solid line) belongs to the refractive 
index of the X-mode and An 2 , (dashed line) to the one of the 
O-mode. The assumed para meters are those of the canonical 
pulsar defined in Sect. 2.2.4. For the O-mode, the magnitude 
of An 2 is much larger and the frequency dependence is much 
stronger than for the X-mode. The difference between the two 
refractive indices clearly decreases with increasing frequency. 



and superposed (leading to weakly polarized parts of the 
profile with often erratic variations in the PPA). 

A spatial separation of polarized modes in a medium is 
caused by birefringence, where both modes propagate in- 
dependently with different indices of refraction (a familiar 
example for such a medium is the calcite crystal). Dif- 
ferent indices of refraction lead to an angular separation 
which may result into an independent observation of the 
two modes. The magnitude of this angular separation is 
given by 

6* = (29) 
dips 

( Melrose 1979| ). Clearly the separation depends on the dif- 
ference between the two indices of refraction. The larger 
the difference is, the better is the chance, that we observe 
the two modes separately. If the difference is small, we 
expect increasing superposition of the two modes, thus 
leading to a depolarization. 

Figure [l3| shows the frequency dependence of the re- 
fractive indices of the two natural wave modes. The devia- 
tions of the indices from unity are displayed: n± = 1+An± 
(see Eq. 18). The -l-mode corresponds to the X-mode (the 
polarization vector lies in the plane spanned by k and B), 
the —mode to the O-mode (polarization vector perpendic- 
ular to k and B). 

Two conclusions can be drawn from the diagram: 1. 
The refractive index of the O-mode is much stronger influ- 
enced by propagation effects than the one of the X-mode. 
This corresponds to the result of Barnard & Arons (1986), 
that the X-mode can propagate through the magneto- 
sphere nearly undisturbed, whereas the O-mode is tight 



to the bending of the field lines. 2. The difference between 
the two refractive indices decreases with increasing fre- 
quency. This may lead to increasing superposition of the 
two natural wave modes (either intrinsically or through 
the limited spatial and temporal resolution of the obser- 
vation). If the phases of the modes are uncorrelated, this 
will cause a depolarization of the observed polarization. 

If this conclusion is correct - depolarization at higher 
frequencies through superposition of OPMs - an addi- 
tional depolarization envelope is superposed on our pre- 
vious discussion. This depolarization could be responsible 
for the low number of pulsars which show the increasing 
degree of circular polarization towards high frequencies. 

The difference between the two indices depends on 
other pulsar parameters as well. Figure |l4] displays the 
frequency dependence of this difference for various values 
for E. n + — Ti- is larger for larger values of E. This leads 
to a depolarization at higher frequencies compared to pul- 
sars with a low E. The peak in the curve for E = 5 • 10 25 
W is caused by a resonance when the wave reaches the 
Dopplcr shifted gyro frequency. It is questionable whether 
the made approximation is valid at this particular point. 
Following Melrose (1979) and Beskin et al. (1993) we 
therefore neglect the resonance and join the asymptotic 
course on both sides of the peak. 

4. Conclusions 

We have used the approximation of the dielectric tensor 
for the low frequency limit to derive the general properties 
of the natural wave modes in the radio pulsar magneto- 
sphere. The approximation of the dielectric tensor in the 
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Frequency [Hz] 

Fig. 14. The difference between the refractive indices of the 
two natural wave modes n+ — n_ versus frequency for dif- 
ferent values of E. The difference is larger for pulsars with a 
high E which should therefore suffer depolarization effects at 
higher frequencies that low _E-pulsars. The peak in the curve 
for E = 5 • 10 25 W is caused by a resonance when the wave 
hits the Doppler shifted gyro frequency and should therefore 
be neglected (see text). 



given limit implies that the plasma is sufficiently cold. 
Following Melrose (1979) this should be true for pulsar 
magnetospheres . 

Electromagnetic waves, which propagate through a 
plasma at an oblique angle to the magnetic field, adopt two 
orthogonal natural wave modes with individual indices of 
refraction. We identify these modes with the familiar or- 
thogonal polarization modes. It was our aim to derive the 
polarimetric properties of these modes and calculate their 
dependence on different parameters under special consid- 
eration of the various angles involved. 

We regard only one natural mode assuming that we see 
only one mode at a time. This assumption is of course not 
always true, but the existence of highly polarized emis- 
sion in single pulses indicates that at least sometimes it 
is true. If both modes are recorded simultaneously, it will 
lead to a depolarization of the observed radiation. Our 
results can therefore be regarded as upper limits for the 
observed degree of polarization. In the case of a highly lin- 
early polarized natural mode we cannot predict the degree 
of linear polarization one will observe. But if the mode has 
only a low level of linear polarization, a low level for the 
observed radiation is predicted. 

The main uncertainty of our calculations is the radius 
of limiting polarization (i?Lp)- However, as we are mainly 
interested in the qualitative behaviour of the polarization, 
we have used a fixed value for i?Lp without restriction of 
generality. For a different value of i? L p the numbers may 
change but the general behaviour remains. This is even 
true if one moves i?LP to its possible extremes: the light 
cylinder radius and just above the emission region. 



Our results are in excellent agreement with the ob- 
servations: 1. The degree of linear polarization decreases 
towards high frequencies (Fig. [t]). The decrease in lin- 
ear polarization implies an increase in circular polariza- 
tion (within this concept). This is in agreement with 
recent observations which show the existence of some 
pulsars which have this peculiar polarimetric property 
( von Hoensbroech ct al. 1998 ). 2. Given a set of specific 
pulsar parameters, a variation of P leads to a variation 
in the spin down luminosity E. Figure || shows that the 
degree of linear polarization of pulsars with a high E de- 
creases at higher frequencies than for pulsars with a lower 
E. Thus a correlation between the polarization at high 
frequencies and E should exist. Such a correlation has 
been indeed observed. 3. Similar to the above point, we 
have varied the period, keeping all other pulsar parameters 
fixed (Fig. [l0|). The degree of linear polarization for long 
period pulsars is expected to decrease at lower frequencies 
than for short periods pulsars. This should give an inverse 
correlation between the degree of linear polarization at 
high frequencies and the period. Again, such an inverse 
correlation has been observed at a frequency of 5 GHz. 

No investigation could be made for the change of the 
polarimetric characteristics along the different field lines, 
which are observed during a pulse. The reason is that 
the qualitative change depends heavily on i?LP- For fu- 
ture work it therefore seems to be important to find an 
estimate for this distance. 

The difference between the refractive indices of the two 
modes decreases with frequency. As this implies a closer 
propagation of the modes, we expect an increasing super- 
position. This would lead to a depolarization of the ob- 
served radiation at high frequencies. Such a depolarization 
could account for low abundance of pulsars which show an 
increasing degree of circular polarization with frequencies. 

We note again that our calculations are based on the 
two assumptions that the background plasma is cold and 
is dominated by particles with one sign of charge. The ob- 
vious advantage of these assumptions - negligible spread 
of the energy distribution and no strong pair production 
- is their physical and mathematical simplicity. The more 
it is astounding that we can qualitatively reproduce the 
complex properties of pulsar radio polarization. In gen- 
eral both assumptions are not widely accepted. Either 
the energy distribution of the plasma is usually taken to 
be extended (Arons & Barnard 1986; Lyutikov 1998 and 
references therein) or the secondarily produced electron 
and positron pairs are expected to dominate the parti- 
cle content of the pulsar magnetosphere (e.g. Sturrock 
1971, Cheng & Ruderman 1977 and Daugherty & Harding 
1982). However the expressiveness of our model deserves 
further considerations. 
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